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5 BACKGROUND OF THE INVENTION 

The present invention pertains to fabrication of integrated circuits and more particularly to 
methods and systems for analyzing the performance and manufacturability of integrated circuits. 

)f Miniaturization and scaling of integrated circuits have resulted in two conflicting trends. On one 

; ^ hand, margins available in high-performance designs are shrinking, increasing the likelihood that 

1 0 ? j inevitable variations during manufacturing will cause performance violations. On the other hand, 

s both the amount of variation as a fraction of the feature sizes, and the sensitivity of the 

^ characteristics of transistors and interconnects to manufacturing variations is increasing. This 

;;f makes it essential to model and simulate accurately and efficiently the impact of manufacturing 
variations on the performance of integrated circuits. 

15 Manufacturing variations include intra-die variations, hereinafter referred to as "local mismatch", 
and inter-die variations, hereinafter referred to as "global variations". Such intra-die and global 
variations are random and systematic variations in the material composition, or in the processing 
steps that are used to fabricate the integrated circuits. Local mismatch, or intra-die variations, are 
differences in the electrical properties of circuit elements that affect components of the circuits 

2 0 that are fabricated on the same die. For example, two identically drawn transistors fabricated 



next to each other in any particular chip will still present different threshold voltage values. 
Global variations, or inter-die variations, affect integrated circuit devices by causing random 
differences in the electrical properties of circuit elements, such as transistors, resistors, 
capacitors, memory cells, wires, etc., that are fabricated on different chips from the same wafer, 
5 on different wafers, or on different batches of wafers. As circuits are designed assuming i) a 
given nominal value of the electrical properties of its elements and ii) that identically drawn 
elements behave exactly in the same way, both global variations and local mismatch affect the 
actual circuit performance and are essential to determine the yield of the product. 

Local mismatch has not been adequately modeled in the past. Accurate matching of the electrical 
10 m properties of active and passive elements is fundamental for functional and parametric 

jf -j performances of analog and mixed-signal integrated performance (DP) blocks, such as operational 
^ amplifiers (OPAMPs), digital to analog (D/A) and analog to digital (A/D) converters, phase 
'If: locked loops (PLLs), etc. Achieving target functional and parametric yield of analog and mixed- 
signal components frequently represents a major bottleneck for the global time-to- volume 
15 performance of complex very large scale integrated (VLSI) systems. 

Previous work in modeling and analysis of matching properties of electronic devices, such as that 
described in the paper entitled "Matching Properties ofMOS Transistors ", Pelgrom M., 
Dunjnmayer, A, and Welters A, IEEE Journal of Solid State Circuits, Vol SC-24, pp 1433- 
1440, October 1989, which paper is incorporated by reference in this detailed description as if 
2 0 fully set forth herein, aims at deriving a suitable model of intra-die metal oxide semiconductor 
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field effect transistor (MOSFET) variance as a function of device size, layout distance and 
orientation. Although these models have become popular, their direct application in the context 
of electrical circuit simulation is difficult for two main reasons. 

First, they model the matching properties of MOSFET "macro" characteristics, such as threshold 
voltage (V Th) or saturation current (I DSAT ) that are only indirectly related to the actual "low-level" 
parameters of most widely used compact SPICE (System Program for Integrated Circuits 
Emphasis) simulation models, such as BSIM3v3 ("MOSFET Modeling andBSIM User Guide, 
Cheng, Y. andHu, C, Kluwer Academic Publishers, Boston, 1999, incorporated by reference as 
if fully set forth herein) or MOS9 ( "Compact Modeling for Analogue Circuit Simulation", 
Velghe, R., et al, IEDM Tech. Digest, pp. 485-488, 1993, incorporated by reference as if fully set 
forth herein). Therefore, a non-trivial inverse modeling process must be applied to extract the 
proper covariance structure of low-level SPICE model parameters corresponding to the available 
matching characterization data for these macro parameters. 

Second, applying a device level mismatch model to the statistical simulation of electronic circuits 
requires the assumption that every matched device is described by a different set of low-level 
device parameters, each associated with a corresponding random variable (RV) ("Statistical 
Modeling of Device Mismatch for Analog Integrated Circuits", Michael, G and Ismail, M, IEEE 
Journal of Solid-State Circuits, Vol 27, No 2, February 1992; "Applying a Submicron Mismatch 
Model to Practical IC Design", Guardiani, G et ah, IEEE CICC Conference, May 1994; 
"Hierarchical Statistical Circuit Characterization of Mixed-Signal Circuits using Behavioral 



Modeling", FeltetaL, IEEE-ACM International Conference on Computer Aided Design, San 
Jose, CA, November 1996, all of which are incorporated by reference as if fully set forth herein). 
The variance of the relevant circuit performance parameters can then be estimated via Monte- 
Carlo analysis. This process requires the generation of a sequence of correlated vectors of 
random numbers, and the evaluation of the circuit performance corresponding to each random 
vector instance either by directly using SPICE or via RSM macro-modeling. The dimensionality 
of the corresponding RV space can be very large when the simultaneous variation of all matched 
n-tuples of devices is considered. 

The problem associated with the large dimensionality of the mismatch simulation task has not yet 
been properly addressed. The G-space approach of Michael and Ismail (cited above), which can 
be proven to be equivalent to the Choleski factorization technique used by Felt et al. (cited 
above), requires j£ (N(mj) - 1) X N(pj) different RVs, where, N( rrij) is the number of matched 
devices of type;, N(pj) is the number of independent process factors used in the model of the; th 
device type, and nd is the number of different devices in the circuit. 

The empirical approach of Guardiani et al (cited above), has an even greater complexity, and can 
be only applied to very simple circuits. Conti {"Parametric Yield Formulation ofMOS ICs 
affected by mismatch effect", Conti, M., IEEE Transactions on Computer Aided Design, vol 18, 
pp. 582 -596, May 1999, incorporated by reference as if fully set forth herein) proposed a method 
based on the experimental characterization of a parametrized auto-correlation function for the 
relevant process parameters described as spatial stochastic processes. The autocorrelation 



function is then used to derive a symbolic formula for the system covariance matrix as a function 
of the layout parameters. Therefore, the complexity of this methodology is also proportional to 
the same number of variables as the G-space approach, however this technique is compatible 
with the statistical simulation methodology described herein, and can be used to replace the 
5 approach of deriving the component correlation matrix using a mismatch model 

SUMMARY OF THE INVENTION 

A method of modeling and simulating the impact of global and local variation on the 
performance of integrated circuits includes the steps of estimating a representation of component 
p mismatch from device performance measurements in a form suitable for circuit simulation; 
1 0 =-1 reducing the complexity of statistical simulation by performing a first level principal component 
{" £ or principal factor decomposition of global variation, including screening; further reducing the 
f «: complexity of statistical simulation by performing a second level principal decomposition 
r ; including screening for each factor retained in the step of reducing the complexity; and 
O performing statistical simulation with the joint representation of global and local variation 
1 5 obtained in the step of further reducing the complexity. 



A method of modeling and simulating the impact of global and local variation on the 
performance of integrated circuits includes the following steps. First, convert, if necessary, the 
device mismatch model at the device performance level into a model suitable for circuit 
simulation, that is, at the SPICE parameter level. Second, if the global variation is provided, 
2 0 perform a first level principal-component or principal-factor decomposition and screening to 
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represent the global variation. If no global variation is provided, perform this decomposition for 
the local variation alone. Third, perform a second level principal-component or principal-factor 
decomposition and screening for each independent factor identified in the second step to 
represent the local variation and reduce its dimensionality. If no global variation is explained in 
5 the initial model before local variation is applied, this step does not apply. Fourth, perform 

statistical circuit simulation and analysis with the combined set of independent factors that result 
from the second and/or third steps to estimate the impact of global and local variation jointly 
and/or separately on the circuit of interest. 

] BRIEF DESCRIPTION OF THE DRAWINGS 

10 Figure 1 is a schematic of a low- voltage opamp used for illustrating the method of the present 
f j invention. 

Figure 2 is a histogram showing the impact of mismatch on the input offset voltage of the opamp 
shown in Figure 1. 

Figure 3 is a block diagram of a digital to analog converter used for illustrating the method of the 
15 present invention on circuits involving matching of more that two components at a time. 
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Figure 4 is a graph showing the efficiency of principal component analysis followed by screening 
in reducing the complexity of statistical simulation including global and local effects. 



In accordance with the present invention, an efficient representation of both inter-die and intra- 
die variation is constructed to analyze the joint impact of these sources of process variation on a 
design. The impact of manufacturing variations on a component is typically represented by 
estimating the distribution of SPICE model parameters for that component. A statistical SPICE 
model that represents both inter-die and intra-die variations has to account for two types of 
correlation; that is, correlation between model parameters and correlation between matched 
components. The correlation between model parameters arises because most commonly used 
SPICE models utilize non-independent parameters. The correlation between different 
components on the same die arises because of intra-die process variation. 

Both the inter-die and intra-die correlation can be represented in a single correlation matrix. 
Given n matched components (C ; „..,C n ) and a SPICE model for each component with m parame- 
ters (P„...P m ), a straight-forward implementation would construct anmxnm correlation matrix: 



DETAILED DESCRIPTION OF THE INVENTION 



1 P 



11,12 
1 




Pl2,ll 



(1) 



PnmM - P 



«m a «(w-l) 



where the rows and columns are Tp p representing parameter j for transistor z. The element p tji 



represents the correlation between Tp } and TjPj. 



This representation is computationally very expensive since the matrix grows as 0(nm x nm). 
Two simplifications are realized by using principal component (PC) decomposition. The first 
simplification is the use of PCS to represent the correlation between the model parameters and 
the second simplification is the use of PCS to represent the correlation between the matched 
devices. 



Principal component representation expresses a set of correlated random variables in terms of a 
set of independent random variables. This is obtained by applying a congruence transformation 
of the form: TTtF 1 = A to the covariance (or correlation) matrix 2, such that T is an orthogonal 
matrix and A is diagonal matrix, its elements k n are the eigenvalues of the covariance matrix. 
An important property of principal component decomposition is that: 



represents the fraction of the variance explained by the first q PCs. By replacing the model 
parameters (p 2 ^p n ) by PCS (£,...,/*) in the Eq. (1) one gets znkxnk size matrix, where typically 
k « m. The reduction in the matrix size stems from the application of a suitable threshold 
filtering algorithm based on Eq. (2). However, further simplifications are possible because using 
PCS instead of model parameters makes the correlation matrix very sparse. Because the PCS are 



Inter-die Principal Component Decomposition 




(2) 
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independent, the correlation between different factors for the same component is zero, i.e., p ilm = 
0 for / * m. Moreover, since the PCS are considered to represent independent sources of 
variations, correlation between different parameters on different transistors is also taken to be 
zero., i.e., p ijkl = 0 for i * k and; # I This leaves only one set of non-zero entries in the 
5 correlation matrix between the same factor for different matched components, i.e., p Ujl . 

Intra-die Principal Component Decomposition 

The correlation matrix is further simplified by using a second level of PC decomposition, replac- 
ing a PC of the original model, e.g., jj, by a linear combination of a set of second-level PCS 
representing the intra-die variation. The result of the second PC decomposition is a unified 
1 0 : ;-i representation that captures both the correlation between the model parameters and the matched 
t components. Moreover, the transformation maintains the correlation between the model 
5 parameters because the RV for each original PC is replaced by an equivalent equation with the 
: : ; same mean and variance. The second level PC transformation adds, at worst, n new PCS, 
t i resulting in a total ofnxk PCS. Since the statistical SPICE model is independent of the 
1 5 application, the number of PCS k is fixed, and usually k « m, where m is the number of model 
parameters. This results in a mismatch simulation method that adds only O(n) new RVs for a 
circuit with n matched components. 

An advantage of the approach of the present invention is that in practice the number of factors 
required is scarcely the n x k factors required in the worst case. This happens for two reasons. 
2 0 First, not all k PCS required to capture the correlation of the model parameters are necessary for 
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representing mismatch. Usually, a much smaller number of variables are required. Second, 
because the amount of mismatch is typically very small, the correlation matrix is, in general, 
characterized by a small number of dominant eigenvalues. Therefore, using the equation that 
gives the fraction of total variance explained by the k-th PCS, the variance of the system can be 
well approximated by using a small number of PCS compared to the worst-case n. 

For example, in the asymptotic case of perfectly matched devices, this method automatically pro- 
duces only one PC (corresponding to the dominant eigenvalue) for each independent process fac- 
tor, resulting in no increase in the number of RVs. In the intermediate situations between the 
worst-case and the asymptotic case, the approach of the present invention provides an 
approximation scheme where a small number of PCS can be selected to accurately approximate 
the correlation between the components. 

Statistical SPICE Models with Mismatch 

Mismatch characterization typically does not produce the correlation between PCS of a statistical 
SPICE model. Usually, the result of mismatch characterization is a set of coefficients for a mis- 
match model of device performance. For example, MOS transistor mismatch characterization 
often results in the coefficients of a Pelgrom-style model for threshold voltage (V^) and satura- 
tion current (I DSAT ) or the transistor gain-factor (£') (see the Pelgrom, et al reference cited above). 
To utilize this information in the representation described above it must be converted into corre- 
lation between PCS of the statistical SPICE model. 
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The procedure for determining the unified statistical representation starts with a statistical SPICE 
model for correlation between model parameters. A subset of the PCS of this model is selected 
for representing component mismatch. The selection is based on two considerations: the device 
characteristics for which mismatch has been characterized, and the weightings of the different 
5 PCS in the equation for each SPICE parameter that impact these device characteristics. For 

example, if mismatch characterization has been performed for threshold- voltage and gain-factor 
of longchannel MOSFETs, then PCS that have the most impact on the VTHO and UO (mobility) 
parameters of the BSIM3v3 model are selected. Once the PCS of the model without mismatch 
have been selected, a numerical optimization is performed to find the correlation between these 
1 0 *i factors in order to obtain the measured mismatch in device characteristics. 

I J Each optimization results in the correlation of the selected PC for one set of device geometry and 
layout parameters. Repeating this procedure for different geometries and layout distances results 
^ in a set of correlation values that are either fit with an interpolating function or represented in a 
r ; : look-up table. This procedure results in a model of PC correlation as a function of device 
1 5 geometry and layout. The step of modeling PC correlation as a function of device geometry and 
layout is performed only once for a particular technology and layout style. 

The method of the percent invention can be used as part of the Circuit Surfer statistical design 
and verification environment. This environment supports many tasks typically required for 
analog and mixed-signal design for manufacturability, such as statistical simulation, sensitivity 
2 0 analysis, response surface modeling, and circuit optimization for manufacturability. Mismatch 
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simulation is implemented in this environment as an annotation of the circuit netlist to specify 
the matched components, and annotations to the statistical SPICE models to include the effect of 
mismatch. The annotated netlist and SPICE models are used to derive a separate statistical 
SPICE model for each component using the two-level PCA. The modified netlist forms the input 
to Circuit Surfer. This implementation enables all the capabilities of statistical design for 
mismatch analyses such as variable screening, response surface modeling and Monte-Carlo using 
mismatch factors, and optimization of a design to reduce its mismatch sensitivity. 

The use of the present invention is illustrated with two applications: a low voltage OP AMP and a 
D/A Converter. These topologies were implemented based on their continued significance in 
low voltage VLSI signal processing applications and their performances sensitivities to device 
mismatch characteristics. 

Low Voltage OPAMP 

The following is a description of an embodiment of the method, in accordance with the present 
invention, for statistical modeling and simulation of the impact of global variations and local 
mismatch on the performance of operational amplifiers (OPAMPs) of the type described. Such 
local mismatch and global variations are typically caused by random and systematic variations in 
the material composition, or in the processing steps that are used to fabricate the integrated 
circuits. Global variations affect integrated circuit devices by causing random differences in the 
electrical properties of the circuit elements, such as transistors, resistors, capacitors, memory 
cells, wires, etc., that are fabricated on different chips from the same wafer, on different wafers, 
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or on different batches of wafers. The local mismatch, or intra-die variations, are differences in 
the electrical properties of circuit elements that affect components of circuits that are fabricated 
on the same die. For example, two identically drawn transistors, fabricated next to each other in 
any particular chip, will still present different threshold voltage values. 

The OP AMP is of the 2-stage rail-to-rail class-AB architecture shown in Figure 1 which employs 
a constant-g m input stage with tail current control. The commonality of complimentary input 
stages in contemporary designs makes this topology ideal for exploring the effects of mismatch 
on offset voltage (Vos) given both NMOS and PMOS mismatch effects. In addition, low voltage 
design examples such as this help exploit the increasing dependence of proper signal resolution 
on low Vos. A Pelgrom model for the mismatch in the threshold voltage and gain factor of long- 
channel transistors is available for this class of technologies. For the matched NMOS pair Ml- 
M2, the model specifies 0(A(Vth)) = 2mFand G(A(£')) = 0.7 %. For the PMOS pair M3-M4, the 
mismatch was: (a(A(Vfft)) = 1.5m V and a(A(£')) = 0.45 %, where Vth represents the MOSFET 
threshold voltage and k' represents the MOSFET gain factor. 

A statistical SPICE model was estimated for the process technology used in this example. The 
model was estimated using measured variance and correlation between key transistor 
performances. This set of performances includes: transistor drive current in the linear and 
saturation regions (Idsat, Idlin) and transistor threshold voltage in the linear and the saturation 
region (Vtsat, Vtlin). The statistical SPICE model is represented in terms of principal factors, 
thus achieves the effect of first-level principal component decomposition. 
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The mismatch model described above is not suitable for direct use in circuit simulation. A 
SPICE level representation of mismatch is derived by selecting SPICE model parameters that 
effect long-channel threshold voltage and the gain factor. For BSIM3v3 SPICE models, these 
parameters are VTHO, TOX, and UO. Numerical optimization is performed to find the 
5 correlation between these factors in order to obtain the measured mismatch in threshold voltage 
and gain factor. 

The statistical model including mismatch was used to simulate the impact of mismatch on key 
O performance measures on this design. This simulation makes use of the method of this invention 
V = to perform the second level principle component decomposition. The impact of mismatch on this 
1 0 design is shown in Table 1. 

Table 1: Impact of Mismatch on OPAMP (10,000 sample Monte- Carlo) 



Performance 


Without Mismatch 


With Mismatch 


Mean 


std. dev 


Mean 


std. dev 


Vos (mV) 


-1.803 


0.09627 


-1.802 


0.7216 



A unified model is derived by applying the method described in the previous sections. The first 
1 5 level PCA results in nine first level PCS (pcal...pca7, pcaLint and pcaWint) for the four 

matched transistors, thus resulting in a 36 x 36 sparsified covariance matrix. A second level 
PCA is performed using this matrix. Three dominant eigenvalues explain over 98% of the total 
variance. In this way, it is possible to account for both inter-die and intra-die variations with the 
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addition of three variables. These three variables are suitably obtained by re-defining pca2, pca4 
and pca5 as linear combinations of two independent components. 

Two sets of Monte-Carlo simulations were performed, one using only the inter-die model and the 
second using the unified model. As expected, this design shows that Vos is extremely sensitive to 
5 mismatch. By quantifying the exact impact of mismatch the method of the present invention 
allows an accurate assessment of manufacturability of this topology. It also shows that 
parametric yield estimates can be overly optimistic for designs sensitive to mismatch if mismatch 
O effects are not accounted for in statistical simulation or worst-case models. Figure 2 shows the 
^ ; distribution of Vos with and without mismatch. The overlay illustrates the 7X increase in the 
1 0 Kl offset voltage standard deviation. Without the ability of efficient and accurate mismatch analysis 
= j the impact on parametric yield due to this increase would be missed. 

; | D/A Converter 

The D/A Converter shown in Figure 3 uses a binary weighted current implementation with eight- 
bit resolution. The output current is summed via current switching and output to a linear 
1 5 current/voltage (I-V) converter. Data converter non-linearity is a measure of the error induced by 
the converter and it is sensitive to the matching between the transistors comprising the ratioed 
current source. 
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Simulation with and without intra-die variation was carried out as described for the OPAMP . 
The procedure was deriving the statistical SPICE model using principal-factors and estimating a 
representation of device mismatch in form of a SPICE model are identical to those for the 
OPAMP. The simulation was carried out using the method of this invention. The mismatch 
5 effects were examined for the binary weighted current source identified in Fig. 3 . Specifically, 

the analog output was tested for differential non-linearity (DNL), defined as the deviation of each 
set of adjacent steps at the analog output from their ideal value (1 LSB). Here, we have defined 
the performance specification to be +/- Vi LSB. The results of statistical simulation with and 
; ; = without intra-die effects are shown in Table 2. 

1 0 h - ] Table 2: Impact of mismatch on D/A Converter (10,000 sample Monte-Carlo) 



Performance 


Without Mismatch 


With Mismatch 


Mean 


std. dev 


Mean 


std. dev 


DNL (LSB) 


-0.176 


0.086 


-0.198 


1.155 


Yield 


100% 


34% 



Intra-die variation causes the standard deviation of DNL to increase by a factor of 13. Such a 
1 5 large increase in the increment causes large changes in the linearity of the D/A Converter. As in 
the previous example, including intra-die effects is vital to get accurate estimation of circuit 
performance variability, but since the circuit's performance for this topology depends upon the 
pairwise matching of eight different transistors, a complete statistical simulation using a standard 
methodology would entail running a Monte-Carlo experiment with large number of correlated 
2 0 RVs. 
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For example, a circuit with eight matched transistors and 7 to 12 independent sources of process 
variation (typical of advanced CMOS processes) requires the generation of a sequence of corre- 
lated random numbers for as many as 96 different RVs. Several thousand SPICE evaluations 
would be required to stabilize the Monte-Carlo results for such a high-dimensional system. The 
5 advantages of complexity reduction using second order principal component decomposition are 
illustrated by the graph depicted in Figure 4. This figure shows the standard deviation of DNL 
versus the number of additional factors employed for mismatch simulation. In fact, although in 
the worst case the present technique would require 18 additional RVs, by applying the screening 
methodology previously described under the heading "Intra-die Principal Component Decompo- 
1 0 \ I sition" above, it is possible to model very accurately the impact of intra-die variability on the 

D/A Converter with only two additional variables. This result makes the simulation of mismatch 
n applicable to an inherently larger set of matched transistors; a requirement as supply levels and 
* feature sizes continually decrease and previously negligible device effects become increasingly 
significant. 
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CLAIMS 



We claim: 



1 . A method for statistical modeling and simulation of the impact of global variation and 
local mismatch on the performance of integrated circuits, comprising the steps of: 

a) estimating a representation of component mismatch from device performance 
measurements in a form suitable for circuit simulation; 

b) reducing the complexity of statistical simulation by performing a first level 
principal component or principal factor decomposition of global variation, 
including screening; 

c) further reducing the complexity of statistical simulation by performing a second 
level principal component decomposition including screening for each factor 
retained in step b to represent local mismatch; and 

d) performing statistical simulation with the joint representation of global variation 
and local mismatch obtained in step c. 



2. A method of modeling and simulating the impact of local mismatch on performance of 
integrated circuits comprising the steps of: 

a) estimating a representation of component mismatch in a form suitable for circuit 
simulation from device performance measurements; 

b) reducing the complexity of statistical simulation by performing principal 
component or principal factor decomposition for local mismatch, including 
screening; and 
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8 c) performing statistical simulation with local mismatch obtained in step b. 

1 3. The method of claim 2 where the principal component or principal factor decomposition 

2 is performed using eigenvalue eigenvector decomposition. 

1 4. The method of claim 1 where the first principal component or principal factor 

2 decomposition is performed using eigenvalue eigenvector decomposition. 

1 _ 5. The method of claim 1 where the second principal component or principal factor 

2 IX decomposition is performed using eigenvalue-eigenvector decomposition. 

1 if J 6. A method for statistical modeling and simulation of the impact of global variation and 

2 local mismatch on the performance of integrated circuits, wherein said method is 

3 «;I integrated in a statistical design and optimization computer-aided design tool to perform 

4 r ; statistical simulation of joint and separate impact of global variation and local mismatch 

5 on performance of integrated circuits and said method comprises the steps of; 

6 a) estimating a representation of component mismatch from device performance 

7 measurements in a form suitable for circuit simulation; 

8 b) reducing the complexity of statistical simulation by performing a first level 

9 principal component or principal factor decomposition of global variation, 
1 0 including screening; 
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11 c) further reducing the complexity of statistical simulation by performing a second 

12 level principal component decomposition including screening for each factor 

13 retained in step b to represent local mismatch; and 

14 d) performing statistical simulation with the joint representation of global variation 

1 5 and local mismatch obtained in step c. 

1 7. A method of modeling and simulating the impact of local mismatch on performance of 

2 integrated circuits, wherein said method is integrated in a statistical design and 

3 r 3 optimization computer-aided design tool to perform statistical simulation of joint and 

4 i "< separate impact of global variation and local mismatch on performance of integrated 

5 k J circuits and said method comprises the steps: 

6 I ] a) estimating a representation of component mismatch in a form suitable for circuit 

7 p simulation from device performance measurements; 

8 r -? b) reducing the complexity of statistical simulation by performing principal 

9 j f component or principal factor decomposition for local mismatch, including 

10 screening; and 

11 c) performing statistical simulation with local mismatch obtained in step b. 

1 8. The method of claim 7 where the principal component or principal factor decomposition 

2 is performed using eigenvalue eigenvector decomposition. 
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The method of claim 6 where the first principal component or principal factor 
decomposition is performed using eigenvalue eigenvector decomposition. 

The method of claim 6 where the second principal component or principal factor 
decomposition is performed using eigenvalue-eigenvector decomposition. 
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ABSTRACT 

A method of modeling and simulating the impact of global variation and local mismatch on the 
performance of integrated circuits includes the step of converting, if necessary, the device 
mismatch model at the device performance level into a model suitable for circuit simulation; for 
example, at the SPICE parameter level. Second, if global variation is provided, performing a 
first level principal-component or principal-factor decomposition and screening to represent the 
global variation. If no global variation is provided, this decomposition is performed for the local 
variation only. Third, performing a second level principal-component or principal-factor 
decomposition and screening for each independent factor identified in the second step to 
represent the local variation and reduce its dimensionality. If no global variation is provided in 
the initial model before local variation is applied, this step does not apply. Fourth, performing 
statistical circuit simulation and analysis with the combined set of independent factors that result 
from the second and/or third steps to estimate the impact of global and local variation jointly 
and/or separately on the circuit of interest. 

PH1Y761963.1 
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Impact of mismatch on offset voltage of an opamp 
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IN THE UNITED STATES 
PATENT AND TRADEMARK OFFICE 

Declaration and Power of Attorney 



As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated below next to my 

name. 

I believe I am the original, first and joint inventor of the subject matter which is 
claimed and for which a patent is sought on the invention entitled EFFICIENT METHOD 
FOR MODELING AND SIMULATION OF THE IMPACT OF LOCAL AND GLOBAL 
VARIATION ON INTERGRATED CIRCUITS the specification of which is attached 
hereto. 

I hereby state that I have reviewed and understand the contents of the above 
identified specification, including the claims, as amended by an amendment, if any, 
specifically referred to in this oath or declaration. 

I acknowledge the duty to disclose all information known to me which is material 
to patentability as defined in Title 37, Code of Federal Regulations, 1 .56. 

I hereby claim foreign priority benefits under Title 35, United States Code, 1 19 of 
any foreign application(s) for patent or inventor's certificate listed below and have also 
identified below any foreign application for patent or inventor's certificate having a filing 
date before that of the application on which priority is claimed: 

None 

I hereby claim the benefit under Title 35, United States Code, 120 of any United 
States application(s) listed below and, insofar as the subject matter of each of the 
claims of this application is not disclosed in the prior United States application in the 
manner provided by the first paragraph of Title 35, United States Code, 1 12, 1 
acknowledge the duty to disclose all information known to me to be material to 
patentability as defined in Title 37, Code of Federal Regulations, 1 .56 which became 
available between the filing date of the prior application and the national or PCT 
international filing date of this application: 

Provisional Application No. 60/166,242 Filed 11/18/99 

I hereby declare that all statements made herein of my own knowledge are true 
and that all statements made on information and belief are believed to be true; and 
further that these statements were made with the knowledge that willful false 
statements and the like so made are punishable by fine or imprisonment, or both, under 
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Section 1001 of Title 18 of the United States Code and that such willful false 
statements may jeopardize the validity of the application or any patent issued thereon. 

I hereby appoint the following attorneys with full power of substitution and 
revocation, to prosecute said application, to make alterations and amendments therein, 
to receive the patent, and to transact all business in the Patent and Trademark Office 
connected therewith: 



Attorney Name(s): William H. Murray Reg. No.: 27,218 

PeterJ.Cronk 32,021 

Robert E. Rosenthal 33,450 

N. Stephan Kinsella 37,657 

Steven E. Koffs 37,163 

Darius C. Gambino 41 ,472 

Lewis F. Gould. Jr. 25,057 

Steohan P. Gribok 29,643 

Samuel W. Apicelli 36,427 

Anthony Colesanti 42,428 

Richard A. Paikoff 34,892 

Gail A. Dalickas 40,979 

Richard T. Redano 32,292 

Joseph A. Powers P47,006 



Telephone calls should be made to William H. Murray 

Phone No.: f215) 979 -1264 
Fax No.: f215^ 979-1020 



All written communications are to be addressed to: 



William H. Murray 
Customer #8933 
DUANE, MORRIS & HECKSCHER LLP 
Customer No. 8933 
One Liberty Place, Suite 4200 
Philadelphia, PA 19103-7396 
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Inventor's signature 




Residence: 7600 Sandia Loop 
Austin, TX 78725 

Citizenship: USA 

Post Office Address: Same as above 

Full name of 4 th inventor: Patrick D. McNamara 
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Residence: 560 N. 6 th Street 
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San Jose, CA 95112 

Citizenship: USA 

Post Office Address: Same as above 
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Residence: 1 240 Klee Court 
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Citizenship: USA 

Post Office Address: Same as above 
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